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Abstract 

We have developed the Smoothed Bandpass Cahbration (SBC) method and the 
best suitable scan pattern to optimize radio spectroscopic observations. Adequate 
spectral smoothing is applied to the spectrum toward OFF-source blank sky adjacent 
to a target source direction for the purpose of bandpass correction. Because the 
smoothing process reduces noise, the integration time for OFF-source scans can be 
reduced keeping the signal-to-noise ratio. Since the smoothing is not applied to ON- 
source scans, the spectral resolution for line features is kept. An optimal smoothing 
window is determined by bandpass flatness evaluated by Spectral Allan Variance 
(SAV). An efficient scan pattern is designed to the OFF-source scans within the 
bandpass stability timescale estimated by Time-based Allan Variance (TAV). We 
have tested the SBC using the digital spectrometer, VESPA, on the VERA Iriki 
station. For the targeted noise level of 5 x 10~^ as a ratio to the system noise, the 
optimal smoothing window was 32 - 60 channels (ch) in the whole bandwidth of 
1024 ch, and the optimal scan pattern was designed as a sequence of 70-s ON + 
10-s OFF scan pairs. The noise level with the SBC was reduced by a factor of 1.74 
compared with the conventional method. The total telescope time to achieve the 
goal with the SBC was 400 s, which was | of 1200 s required by the conventional 
way. Improvement in telescope time efficiency with the SBC was calculated as 3x, 
2x and 1.3 x for single-beam, dual-beam, and on-the-fly (OTF) scans, respectively. 
The SBC works to optimize scan patterns for observations from now, and also works 
to improve signal-to-noise ratios of archival data if ON- and OFF-source spectra are 
individually recorded, though the efficiency depends on the spectral stability of the 
receiving system. 

Key words: telescope — technics: calibration — techniques: spectroscopy — 
instrumentation: detectors — methods: statistical 



1. Introduction 



Spectral observations in radio astronomy are aimed to detect emission or absorption 
line features that bring us rich information such as intensity, velocity, line width, etc. To 
acquire a spectrum, signals received by a radio telescope are processed in a spectrometer. The 
performance of a spectral observation is characterized by the spectral resolution, the bandwidth, 
the sensitivity, and the stability. 

A certain bandpass calibration procedure is needed to obtain a desired performance 
of spectral observations. An acquired spectrum is a summation of unwanted noise, Tgys, and 
the signal from the target source that is denoted as the antenna temperature, Ta(z/). The 
spectral shape is affected by the bandpass response, H{u), that is the transfer function of the 
receiving system as a function of frequency, u. The acquired ON-source spectrum will be then 
^on(i^) = H{h')(Tsys + Tg.ii')). Bandpass calibration is necessary to estimate Ta(z/) from the 
acquired spectrum, Ton(i^)- 

In most of radio observations, the acquired spectrum is dominated by the system noise, 
i.e., Tgys ^ Ta(i^). To relieve the source spectrum from the system noise, a spectrum of OFF- 
source blank sky adjacent to the target source is subtracted from the ON-source spectrum. 
This position-switching scan yields the source spectrum as 

^ _ ToN(t^)-ToFF(t^) ( . 

Jal^/j- ^^^^ , llj 

where Tqff is the spectrum taken from the OFF-source scan. In conventional position-switching 
observations, the bandpass Hlu) is also acquired from the OFF-source scans because the blank 

sky emits a featureless flat spectrum. For H{iy) = — — , equation 1 will be 

-'sys 

The sensitivity of spectal observations is evaluated by a SD (standard deviation), a, of 
the acquired spectrum in line-free channels. It is derived from equation 2 as 

\T^{iy) + Tsys J \Ton{i^)J \Toff{i^>/ 
where cton and ctoff are SDs of the ON- and OFF-source spectra, respectively. For weak 
sources for which To^{i^) — Toff(i^), a is given by a root of sum squared (RSS) of cton and 
(JoFF- When cton and ctoff are dominated by thermal noise, according to Wilson et al. (2009), 
they are expected to relate to the integration time, toN and toFF, the spectral resolution, u^cs, 
and Tgys as followed, 

T -\-T T 
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The total telescope time is the sum of toN, ^off, and other overhead time such as setup, system 
noise measurement, gaps between scans, etc. Optimized scan pattern is designed to maximize 
the sensitivity, i.e., to minimize a, within the minimal telescope time. Conventionally toN and 
toFF are equally allocated in position-switching observations to have cton — o"off- This manner 
makes the total telescope time longer than twice of toN and lets a ~ -\/2croN- 

If there is an efficient way to reduce ctoff less than o"on, we can save the total telescope 
time and get better sensitivity. For an ideal receiving system, whose bandpass response were 
known a priori and very stable, OFF-source observations would be unnecessary. Modern sys- 
tems offer better spectral stability in the trend in replacing analog devices with a high-speed 
digital sampling system in an upper stream of signal transfer. The bandpass response of a 
digital system, including a digital filter and a digital spectrometer, is determined by a fixed 
algorithm and is not affected by environmental variations. Utility of digital devices brings us a 
possibility to reduce the integration time and the SD of OFF-source scans. 

In this paper we propose the Smoothed Bandpass Calibration (SBC) method where the 
OFF-source spectrum is smoothed across the bandwidth to reduce the integration time of the 
OFF-source scans and the SD of the spectrum. Section 2 describes the concept of the SBC 
method. Section 3 reports test experiments that verify the SBC and that quest for the optimal 
parameters of a smoothing window and a scan pattern. In section 4 we discuss the behavior of 
the spectral stability and estimate sensitivity improvements expected in some scan patterns. 

2. Method 

We propose relevant spectral smoothing along frequency for OFF-source scans to reduce 
required integration time that achieves a desired SD. Before the spectral smoothing, a bandpass 
correction is employed using a template bandpass response that is presumed to be relevantly 
static and taken by long-time integration toward blank sky. The bandpass-corrected OFF- 
source spectrum will be almost fiat with slight warp due to fiuctuations of the receiving system 
and smeared by thermal noise. In the following subsections we formularize the spectral behavior 
and establish a strategy to reduce spectral unevenness. 

2.1. Description and Assumption 

As described in section 1, the purpose of OFF-source scans is to subtract Tgys and 
calibrate the bandpass, H{u), which is a frequency- dependent system response. Since the 
bandpass can be time-variable, we denote it as H{v.,t) = Hq{v) + Hi{i',t), where Hq{v) is a 
static template bandpass and Hi{v^t) is a fiuctuation from Hq{v). The observed OFF-source 
spectrum is composed of not only H{v.,t) but also thermal noise, e(z/, t), which is assumed to 
be attributed to a stochastic random process following the normal distribution with the mean, 
/i = 0, and the variance, cr'^oisc — ' where tintog is the integration time. Thus the observed 

^res^integ 

OFF-source spectrum will be ToffI^^,^) = {.Hq{u) + iJi(z/,t))Tsys(l + e(z/,t)), and its expectation 
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will be (Toff('^,^)) = {Ho{iy) + Hi{iy,t))Tsys- Sufficiently long integration time is required for 
conventional OFF-source scans to reduce dnoise- 

We assume that Hi{v,t) is relatively smaller than Hq^u), almost fiat and smooth along 
the frequency, with a typical frequency scale, z/var of variation. Relevant spectral smoothing for 

OFF-source spectra efficiently reduces cinoisc by a factor of W— ^ and keeps the expectation of 

OFF-source spectrum. Since the spectral smoothing will not be applied for ON-source spectra, 
desired spectral resolution will be kept. 

2.2. Allan Variance Decomposition 

The flatness of Hi{v.,t) is evaluated by the spectral Allan variance (SAV) analysis. The 
SAV in defined as 

[H{v + Az/) - 2H{y) + H{v - Az/)]2 ' 
2Ai^2 

For a small frequency span (Az/ <C z/^ar) the SAV is dominated by thermal noise and follows the 
characteristics of cry(Az/) oc Az^~^. At a larger frequency span (Az/ ^ z/var)? Hiiy.t) exceeds the 
thermal noise and the SAV shows different characteristics from that of the thermal noise. SAV 
analysis allows us to decompose ifi(z/,t) from the thermal noise. 

An example of SAV for an OFF-source spectrum is shown in figure 1, as a function of the 
channel spacing, z/gp^. The SAV in the 1-min integrated raw spectrum (a) is decomposed into 
two components; (1) thermal noise, e(z/,t), that dominates in the range of ^<5and whose SAV 
has a dependence, Az/~^ and (2) bandpass characteristics, H{i>,t), that dominates in the range of 
^ > 5. When components (1) and (2) show approximately symmetric power-law characteristics, 
the bottom of the SAV appears at the frequency span where those two components are almost 
equivalent (see results in subsectino 3.2 and discussion in subsection 4.1). After bandpass 
correction using a template spectrum, HqIu), component (2) is considerably eliminated as (b). 
The rest fluctuation, Hi(iy,t), dominates the SAV in ^ > 120. As the thermal noise will be 
reduced via time integration in (c), the bottom of the SAV shifts to smaller Az/ (e.g. ^ ~ 40 
for 546- min integration). Relevant spectral smoothing with the smoothing window given by 
the bottom in SAV, as described in the next subsection, will yield efficient reduction of the 
thermal noise keeping the bandpass response as shown in (d). 

2.3. Strategy for Adequate Smoothing and Scans 

The relevant smoothing window and ON-OFF scan pattern should be determined to 
minimize the SD in the OFF-source-subtracted spectrum within the given telescope time, or 
to minimize the total telescope time that achieves the given SD. The best smoothing window, 
A'sw = — is determined by the bottom of the SAV. After we apply spectral smoothing to the 



As described in section 3, we used an FX spectrometer whose spectral resolution function is a squared sine 
function and i/gp = v^^s by 2048-point FFT to produce 1024 ch. 
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OFF-source spectra, their variance is expected to be reduced to ^iB^. 

''sw 

Because of inequality in variances of ON- and OFF-source spectra, we need to re-design 
the best duty cycles of toN and toFF under the condition where total telescope time is constant. 
Let an ON-to-OFF-source integration time ratio set to 1 — x : x in each of the duty cycle and 
the total telescope time of ttot^- Substituting toN = (1 —x)ttot, ^off = ^^^^tot, and assuming that 
^ Tsys — Ton — Toff in equations 3 and 4, we have 

a \ 1 



+ ^ • (6) 



The variance is minimized to 

/ \2 
/ CT rn\r\ \ -L 



^-1!. (7) 



\ -^sys / ^res^tot ^sw 

when we have x = , When A'^.^ = 45 for instance, the best ratio is 0.13. 

The interval between adjacent OFF-source scans should be shorter than the timescale 
of spectral stability in the receiving system. The spectral stability is evaluated by the time- 
based Allan Variance (TAV). As we'll see in figure 3, a typical spectral shape of Hi{u,t) show 
an S-shaped feature with a peak and a bottom. Since we focus on stability of the bandpass 
shape, the TAV is calculated by the difference between levels of the peak and the bottom, i.e. 
AHi{t) = Hi{upcak,t) - Hi{ukottom,t). Thus, the TAV is derived as 

/ [AH, it + r)-2AH^it)+AH^{t-T)]' \ 
^y(^) = [ ^2 /' (8) 

where r is the time lag. The stability timescale is determined by the time lag where the TAV 
follows cry(r) oc r~^. 

In summary, our strategy to determine the optimal spectral smoothing window and the 
scan pattern will be: 

1. Acquire the template bandpass response, Hq^u), and its dependence on a frequency band 
by sufficiently long integration before programed observations (e.g. during seasonal main- 
tenance). Relevant spectral smoothing can be applied to the template bandpass to reduce 
the thermal noise in it. 

2. Set the target SD for the observation and estimate approximate integration time. 

3. Acquire test scan data to take SAV and TAV to decide at the bottom of the SAV. 

4. Derive the ON-OFF duty cycle of 1 — x : a; basing on Ns^/. 

5. Estimate the timescale of stability and determine the OFF-source-scan interval. 

6. Design the optimal ON-OFF scan pattern basing on the time stability of the bandpass 
and the ON-OFF duty cycle. 



Here, we omit overhead and scan gaps for the simphcity. 
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7. Apply bandpass calibration of Hq{u) for observed ON- and OFF-source spectra before 
integration. Employ spectral smoothing for the OFF-source spectra with the smoothing 
window of 

8. Apply the bandpass calibration of Hi for ON-source spectra using the adjacent smoothed 
OFF-source spectra. Time-integrate the bandpass-calibrated ON-source spectra. 

9. Apply baseline subtraction, if necessary, to the OFF-source-subtracted spectrum and get 
the final result. 

3. Tests and Results 

Field tests have been carried out for the SBC to investigate how does it work efficiently. 
We used the VERA (Kobayashi et al. 2003) Iriki 20-m antenna in the single-dish mode. The 
antenna pointed to the zenith throughout our tests. The 22-GHz HEMT receiver was tuned 
at 22.235 GHz with the bandwidth of 512 MHz. The received signal was downconverted with 
the first LO of 16.8 GHz, and the second LO of 5179 MHz before digitized at 1024 Msps 2-bit 
quantization by the digital sampler ADS-1000 (Nakajima et al. 2001). Digital filtering was 
applied using the digital filter unit (Iguchi et al. 2005) to split the signal into sixteen 16-MHz 
streams where we used only the first stream. Spectroscopy was employed using the software 
spectrometer, VESPA (VEra SPectrum Analyzer; Beppu et al. 2011), that produces 1024 ch 
by 2048-point FFT. This produces a squared sine spectral resolution function whose first null 
appears at the adjacent ch, i.e., z/gp = z^res- 

The template bandpass, Hq{v), was obtained on 2010 Apr. 15 for 230- min integration. 
A 10-hour observation was conducted on 2010 Oct. 5, half-a-year later than acquiring the 
template bandpass to include seasonal variation. Pseudo ON-OFF scan pattern was produced 
from the 10-hour continuous observation; OFF-source scans was pilfered from the continuous 
zenith observations. Data quality was checked via monitoring the total power of every spectrum. 
The first 50-min data were not used because of unexpected operations of hot-load insertion. 
We also flagged out 3-sec data because of obviously irregular spectra probably due to data 
transmission error. Finally, we obtained 32768-sec continuous spectra. For calculation of SAV, 
we did not use the first spectral channel that includes direct current (DC) component affected 
by voltage bias offset at the digital sampler. 

In this test, we set the target SD of the antenna temperature of 50 mK, which is the 
standard of regular monitoring program in the VERA Iriki single-dish observations. This 
corresponds to tS— = 5 x 10^"^ for T^ys = 100 K. 

sys 

We employed 3rd-order B-spline smoothing for OFF-source spectra in step 7 described in 
subsection 2.3. Node intervals are set to be equal to the smoothing window, Nsy^. The statistical 
package R was used for data analyses of integration, bandpass calibration, AV calculation, and 
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B-spline smoothing. All of source codes for these procedures are presented in the GitHub^. 

3.1. Bandpass Spectra 

The template bandpass, HqIjj) is shown in figure 2. The SD of its random noise was 
1.97 X 10~^, which was estimated by Hq{i') / Hq{v) — 1 where Hq{v) was a B-spline-smoothed 
spectrum. We used Hq^u) for the following tests, though the SD was one order of magnitude 
smaller than the target noise level, to quest the ultimate performance. 

Bandpass-corrected spectra, Q^^^ ' _ shown in figure 3. Here, Tqff is the 

mean value of the OFF-source spectra that represents Tsyg. Thus the bandpass-corrected spec- 
trum indicates - ^ (1 + e(i/,t)) +e(i/,t). The solid line in figure 3 was computed by the 3rd- 

order B-spline smoothing with the node intervals, Ns^ = 45 and represent - ^ ' } . Departures 

from the solid line were dominated by the thermal noise, [ IH -^7— t" I ^i^^t) — ^i^^t)- The 

smoothed spectrum was time variable. Its time variability is described in subsection 3.3. 

3.2. Spectral Allan Variance 

The SAVs of bandpass-corrected spectra defined as equation 5 are shown in figure 4. At 
every frequency channel separation, the SAV decreased as was integrated for longer time. While 
the SAV was decreasing as a function of frequency channel separation at shorter integration 
than 16 s, it appeared a local bottom and a top for longer integration than 16 s. They appeared 
at < 200 ch and ~ 259 ch, respectively. Figure 5 shows the SAV of a 64-min-integrated spectrum 
as an example. The profile was composed of two power-law components. For narrower channel 
separations than the bottom, the SAVs was dominated by a power-law profile with the index of 
— 2.003±0.001. The index above the bottom was 1.29±0.02. The bottom of the SAV appeared 
near the frequency span where the two power-law components were equivalent. 

At longer time integration, the bottom shifted towards narrower spectral separation in 
longer integration time, while the top stayed at almost the same separation. These behavior is 
summarized in table 1. 

3.3. Time Stability 

Figure 6 displays variability of - ^ ' } . The smoothed spectrum shows time variation in 

terms of not only the power level but also its shape such as slope, curvature, and local bumps. 
Variation of the bandpass was evaluated by the standard deviation of the smoothed spectrum, 

(Ttiiy), defined as cr^lu) = I ( -^^ ' ^ ] \ — ( r}/^\ / ■ Here, the expectation is taken by time 
average. 

Figure 7 shows the standard deviation of the smoothed spectrum. Two clear peaks 



https:/ /github.com/kanicnoseiji/BPsmooth 



appeared at 71 ch and 920 ch, near the shoulders of the bandpass, Ho{u). The standard 
deviation of ~ 2% was significantly greater than the target accuracy of = 5 x 10~^. Hence, 
it is necessary to calibrate the bandpass using the smoothed OFF-source spectrum within the 

timescale while / } is stable enough. 

The timescale of stability was evaluated using the TAV defined in equation 8. Since the 
top and bottom channels of the smoothed bandpass were also time variable, we instead used 
the twin peak channels in o"f(z/), i.e. we substituted AHi{t) = ifi(71ch,t) — ifi(921ch,t). The 
TAV is shown in figure 8. The first bottom of 30-sec-integrated spectrum appeared at 300 s. 
The power-law index of the TAV was ~ —2 for r < 60 s and was —1.7 at r = 80, which we 
regulate as the spectral stability timescale in following processes. 

3.4- Scan Patterns and Performance 

Basing on above results we set the new scan pattern and compared with the conventional 
ON-OFF scans without spectral smoothing. The conventional pattern (case 1) consists of pairs 
of 30-s ON and 30-s OFF scans. The new scan pattern (case 2) was a set of 70-s ON and 10- 
s OFF with the SBC. The pattern was designed as the cycle period to be shorter than the 
spectral stability timescale of 80 s and the optimal ON-to-OFF integration time ratio described 
in equation 6 with A^sw = 45. 

To reach the targeted noise level of cr/Tsys = 5 x 10~^, the conventional pattern required 
20 sets (1200 s) of the total telescope time. Its resultant spectrum is shown in figure 9 (Al). 
The case-2 spectrum with the telescope time of 1200 s (15 sets of 80-s scans) is shown in figure 
9 (CI). The SD was (2.8to;2) ^ 10~^ (median, minimum, and maximum of 27 samples), ^jjg 
that of the conventional scan pattern. To achieve the targeted noise level, case 2 required 5 
sets (400 s). Its resultant spectrum is shown in figure 9 (Bl). 

Since the resultant spectrum appeared wiggled features other than random noise, we 
attempted baseline fitting and subtraction using the 3rd-order B-spline function with the node 
interval of 45 ch. After the baseline subtraction, the SDs of cases 1, 2 (1200 s), and 2 (400 s) 
became (4.7lo;2) x 10""^, {2.6lio^) ^ 10"^, and (4.5lo;3) ^ 10"^, respectively, and shown in figures 
9 (A2, C2, and B2). 

The performance of SD as functions of total telescope time is summarized in table 2 and 
figure 10. At any total telescope time, the SBC exhibited better performance than conventional 
scans. While the SD of the conventional scans followed a oc t^;^^ dependence, as was expected, 
the SBC resulted in a shallower slope with the power index of —0.46. 

3.5. Smoothing Window 

The optimal spectral smoothing window, Au, was determined by the bottom in the SAV 
listed in table 1. For the timescale longer than 16 min, the bottom appeared at Au/usp ~ 32 — 60 
ch. We chose the window of 45 ch for time stability tests and spectral performance tests in 
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previous subsections. We also tested the spectral performance with various windows of 2, 3, 4, 
. . ., 191 ch for the scan pattern of case 2 with the total telescope time of 400 and 1200 s. 

Figure 11 shows the results. The median values of SDs in the resultant spectrum before 
baseline subtraction recorded a minimum at the 64-ch window for both telescope times of 400 
s and 1200 s. After baseline subtraction, the SD became significantly lower and the minimum 
appeared a fiat bed for wider window than 45 ch. In all cases the bottom was smaller than the 
targeted noise level. 

4. Discussion 

4.1. Bandpass Flatness 

As is shown in figure 1, the SAV across bandwidth is composed of the thermal noise 
and fluctuation of bandpass characteristics. The crossover point, where the thermal noise and 
the bandpass fluctuation is equivalent, appears at the bottom of the SAV. Relevant smoothing 
window, Au, should be set around the crossover point. For longer integration the bottom of SAV 
shifts toward narrower channel separation. The regression of the crossover point indicates that 
the channel separation of the bottom relates to the integration time as Ai^bottom oc t^ej^^^ '^^^ 
for tintcg < 128 min. 

This behavior is explained as following consideration. The thermal noise component 
is stochastic which does not depend on time and frequency. Thus the SAV of thermal noise 
component oc Au^'^t^^^. The bandpass fluctuation SAV oc Az/^'^^''=° °^. This power-law index is 
slightly steeper than SAV oc Au which is expected for a random-walk process along frequency. 
Since the bandpass characteristics is frozen during short timescale, the SAV of bandpass fluc- 
tuation is independent of tinteg- Therefore, the crossover point appear under the condition of 
^^"^^intcg — «Az/^'^^ where a is a proportional coefficient. Thus, the SAV bottom will appear at 

Az4,ottom octjnteg®. The behavior of the test observation, which indicates Ai^bottom t^iS.cJ^'^^''^^^ ■, 
is consistent with the theoretical consideration. 

Although a different combination of a receiving system and a spectrometer other than the 
VERA system may have different properties of bandpass flatness, measuring SAVs at various 
integration time allows us to decompose it into thermal-noise and random-walk components. 

4-2. Expected Time Reduction 

As shown in subsection 3.4, the SBC (case 2) allows us to shorten the total telscope time 
to I to achieve the targeted noise level, compared with the conventional method (case 1). In 
other words, the SBC reduces the SD by a factor of \/3 within the same total telescope time. 
These results demonstrate that the telescope time efficiency can be tripled by the SBC method 
with the optimal scan pattern. 

Since our analysis did not take scan gaps time into consideration, the efficiency can 
be different in realistic observations. As shown in figure 9, systematic undulation remains in 
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the case-2 spectrum before baseline subtraction. The undulation can be ascribed to OFF-scan 
time interval of 70 s that may be longer than the stability timescale (see figure 8). OFF-scan 
intervals should be shorter to reduce the undulation, however, this modification increases the 
time loss in scan gaps. Optimal scan pattern should be tested in detail under realistic conditions 
of scan gaps and time stability. 

As shown in figure 10, the slope of SD reduction by time integration for the SBC was 
shallower than that for conventional scans. The shallower slope indicates that the SBC under- 
performs at a long integration time. This underperformance can be caused by time instability 
of the bandpass. The power index of the TAV of —1.7 at r = 80 s, with a small departure from 
—2, indicates that bandpass fluctuation remains with the level of Oy ~ 10~^ (see figure 8). 
This level of fluctuation is not tangible in the conventional scans even at the longest integration 
= 9.3 X 10~^) but can affect in the SBC {y~ = 6.3 x 10~^). The bandpass stability can be 
the principal component to determine the sensitivity at a longer integration time then 32400 s. 
To pursue the sensitivity at the longer integration, shorter scan pattern should be considered 
in spite of time loss in scan gaps. 

4-3. Advantages in a Dual Beam System 

Dual beam systems are equipped in some radio telescopes such as the GBT 100-m (Jewell 
2002), the Nobeyama 45-m (Kuno et al. 2010), and the VERA (Kobayashi et al. 2003), to double 
the time efficiency by interchanging ON- and OFF-source scans between two beams. The duty 
cycle for ON and OFF scans in dual beam systems should be 1:1 (or x = 0.5 in equation 7), 
if performance of two receivers is uniform, and different from the optimal duty cycle for the 
SBC. In this case the variance, a^, will be reduced by a factor of - i 1 + — — ) when we employ 
the SBC. Taking Ng^ = 45 for instance, the variance is reduced by the factor of 0.51 and the 
telescope time efficiency is almost doubled, 3.9 times better than conventional ON-OFF scans 
with a single beam system, or 1.3 times better than the SBC with a single beam system. 

Consideration in this subsection implies that the SBC for single-beam with conventional 
ON-OFF scans (1:1 duty cycles) also doubles SNR. Re-analysis with the SBC for previous 
observations can offer an opportunity to enhance SNR, if ON- and OFF-source spectra are 
separately recorded, and the efficiency depends on the spectral stability of the system. 

4.4- Advantages in OTF Scans 

ON-the-fiy (OTF) mapping shares the OFF-source pointing among a continuous scan 
pass (Mangum et al. 2007). Sawada et al. (2008) argued that the noise level achieved in the 
unit observing time of the OTF map is written as 




10 



^OH + ^!^EQ I ^row/cal, (9) 



vscan 



and derived the optimal OFF-source integration time to map a rectangular area with the row 

span of h was given toFF = ]j (4can + ^oh) \J iVscan , where tscan = h/ Vscan is the ON-source 

scan time for a row by the scan speed of Wscan, i^oH is the overhead time between row scans, r] is 
the efficiency determined by the gridding convolution function, d is the grid spacing, and Nf^^ 
is the number of row scans between OFF-source integrations. In the practical case of 77 = 4.3, 
h = 600", Al = 5", d = 7".5, ton = 25 s, and N^^^ = 1, according to Sawada et al. (2008), the 
ON-source scan time of tgcan = 20, 40, and 60 s yields t^p^^^ = 7, 12, and 17 s, respectively. 
The optimization is modified when the SBC is applied. Equation 9 is modified as 



OFF 
sW-cell , 



iVrow/cal, (10) 



vscan 



when the smoothing window is A^^g^ ch. The optimal OFF-source integration will be toFF = 
\/(^scan + ^oh) ^ '^T" \I ^scan . For comparisou to Sawada et al. (2008), tscan = 20, 40, and 60 s 

V ^SW<1 

yields ^qf™'^^ = 1, 1.8, and 2.5 s when we take A^sw = 45, and the noise level decreases by factors 
of 0.886, 0.874, and 0.856, respectively. Thus, the total telescope time that achieves the same 
SD will be reduced by a factor of ~ 1.3 — 1.4. 

The effect of the SBC is less than simple ON-OFF scans, nevertheless, the SBC offers 
somewhat better results. Consideration in this subsection does not involve bandpass fluctuation. 
More realistic optimization and noise level should be estimated basing on time-stability of 
receiving systems. 

4-5. Vulnerability against Spurious Signals 

SBC yields weakness against unwanted spurious emissions such as RFI, contamination 
of sampling clocks and reference signals in a receiving system, artificial pattern caused by 
numerical errors of spectral calculation, etc. These spurious signals usually show spiky features 
in both ON- and OFF-scan spectra. They are suppressed in the smoothed OFF-scan spectra 
while kept in the unsmoothed ON-scan spectra. Thus they will appear in the final results, 
though they would be subtracted by the conventional ON - OFF scans if they were stable. 

In our test observations, spurious emissions appeared at 704 ch and 832 ch (see figure 
9) that correspond to 11/16 and 13/16 of the bandwidth. Their line width is narrower than 
the spectral resolution. Digital noise such as harmonics of the sampling clock can cause such 
spurious emissions whose frequencies are simple fractions of the bandwidth. These spurious 
frequencies should be masked in spectral data reduction. Observers need to pay deeper attention 
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to spurious when they apply SBC. 
5. Summary 

It is presented that the proposed SBC method offers significantly better performance 
than conventional ON-OFF scans, especially when a stable spectrometer such as digital system 
was equipped. Our tests showed that the total telescope time was reduced to | to attain 
the same SNR, or xl.7 better SNR was obtained in the same telescope time for a single- 
pointing ON-OFF scans. The SBC method was also efficient for dual-beam systems, ON-the-fly 
mapping. 

The optimal analysis scheme of the SBC is presented in this paper. To apply this 
method, radio observatories should offer 

• Long-time integrated bandpass to obtain good-enough Hq{u). This can be done in off 
seasons. 

• SAV to show flatness of Hi{u,t) to allow observers to compute the optimal smoothing 
window. 

• TAV to estimate stability timescale of Hi{v,t) to design the optimal scan pattern and total 
telescope time required. 

• Spectral outputs of ON- and OFF-source, separately. Automated ON-OFF subtraction 
makes it impossible to OFF-source spectral smoothing in post-observation analysis. 

Finally, re-analysis of previous spectral data with the SBC can also offer an opportunity 
to enhance SNR, if ON- and OFF-source spectra are separately recorded. 

The VERA Iriki 20-m telescope is a part of the VERA VLBI network which is operated 
by the National Astronomical Observatory of Japan and staffs and students of the Kagoshima 
University. 
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Table 1. Bottom and top appeared in the spectral AV 



Integ. 


bottom 


bottom SAV 


top 


top SAV 


(min) 


(ch) 


X 10-10 


(ch) 


X 10-10 


1 


128 


5.3 


255 


7.2 


2 


103 


4.0 


257 


7.0 


4 


86 


3.1 


259 


6.8 


8 


71 


2.4 


259 


6.7 


16 


60 


1.9 


259 


6.7 


32 


49 


1.4 


259 


6.5 


64 


42 


1.1 


259 


6.5 


128 


39 


0.9 


259 


6.4 


256 


34 


0.7 


259 


6.3 


512 


32 


0.6 


259 


6.3 



Table 2. Spectral performance as functions of total telescope time 



ttot <y/T,y, (xio-i) 



(min) case la case lb case 2a case 2b 



4 


10 


0+3.7 
■^-0.8 


10.8^0;? 


6 


1+0.5 
•-'^-0.4 


5, 


7+0.5 
■ ' -0.4 


8 


7. 


7+2.2 
■ ' -0.5 


7 c,+0.4 


4, 


O+0.4 
•"^-0.3 


4, 


1+0.3 
■-^-0.2 


16 


5. 


4+1.1 
■^-0.4 


c O+0.3 
^■'5-0.3 


3 


1+0.3 
■-■^-0.2 


2, 


q+0.2 
■^-0.1 


32 


3. 


nn+0.28 

■^^-0.25 


O 7r+0.16 

■J- ' ^^-0.20 


2 


oc+O.U 
■^^^-0.08 


2 


07+011 

-0.10 


64 


2. 


7R+0.03 
■ '"-0.14 


2 66+0 04 

^■""-0.10 


1, 


Rq+0.05 
•"«J-0.11 


1, 


49+0.06 

■^^-0.05 


128 


1. 


Qo+0.05 
■^'-'-0.03 


1 aa+0-006 
-l^-^^-o.oie 


1, 


1Q+0.05 
•^^-0.04 


1, 


07+0.02 
■"'-0.01 


256 


1. 


qo+0.01 
■"J'J-0.01 


-1 9Q+0.002 
^•^^-0.003 


0, 


oq+0.01 


0, 


70+0.02 
• '°-0.02 


512 


0. 


,93 


0.91 


0, 


.70 


0, 


.59 



Cases with 'a' and 'b' stand tor before and after baseline subtrac- 
tion, respectively. 
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Wilson, T. L., Rohlfs, K., Hiittemeister, S. 2009, Tools of Radio Astronomy (Berlin: Springer), p. 69 
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1 2 5 10 20 50 100 200 500 

Channel Separation [ch] 

Fig. 1. Spectral Allan Variance (SAV) across the 16-MHz bandwidth with 1024 spectral channels. The 
data conies from the field test described in section 3. Four SAV profiles are derived from (a) 1-min inte- 
grated raw spectrum before bandpass correction, (6) 1-min integrated spectrum after bandpass correction, 
(c) 546-min integrated spectrum after bandpass correction, and (d) 546-min integrated spectrum after 
bandpass correction and spline smoothing are applied. 
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Fig. 2. Template bandpass spectrum, Hq{v), obtained by 230-min integration on 2010 Apr. 15. The DC 
bias (spectrum at ch= 1) is flagged out. The amplitude is normalized as the mean becomes unity. 



S ^ 



400 600 
Frequency [ch] 



Fig. 3. Bandpass-corrected spectrum, 



for 10-s average. The solid line 



shows the smoothed spectrum, / } with the smoothing window of 45 ch. 
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1 2 5 10 20 50 1D0 200 500 

Channel Separation [ch| 



Fig. 4. SAVs across the bandwidth under thue integration of 1, 2, 4, 32768 s. 




1 2 5 10 20 50 100 200 

Channel Separation \ch] 



Fig. 5. SAV at 64-niin integration. The solid line indicates the best fit with two power-law components; 
a^(u) = aiv°'' +a2iy°'^. The best-fit power indices are ai = -2.003 ±0.001 and a2 = 1.29 ±0.02. 
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Fig. 6. Time variation of the bandpass-corrected speetrum, 7^~r- Each spee- 

trum was time-averaged for 30 s and spectral-smoothed with a 45-ch B-sphne window. 




Fig. 7. Standard deviation (SD) 
spectrum, crf(i^) = 



of 




Frequench [ch] 

each spectral channel in the bandpass-corrected 
The spectrum was time-aver- 



aged for 30 s and spectral-smoothed with a 45-ch B-spline window before calculation 
of the SD. Two peaks in the SD appeared at the frequency channels of 71 and 920. 
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100 200 500 
Time Lag [sec] 

Fig. 8. Time-based Allan variance (TAV), <Jy{T), of spectral difference between 71 eh and 920 eh as 
a function of time lag as defined in equation 8. Sequences of filled circles and open squares rep- 
resent 10-s and 30-s time-integrated spectra, respectively. The local TAV minima appeared at 370 
s and 300 s for each. Significant departure from the power-law index of —2 arose at r ^ 60 s. 
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Before baseline subtraction 



After baseline subtraction 




Frequency [ch] Frequency [ch] 

Fig. 9. Resultant spectra after bandpass correction. Left and right panels stand for before and after base- 
line subtraction performed after bandpass correction is applied. Al and A2 : conventional ON-OFF scans 
(case 1) of the total telescope time of 1200 s. Bl and B2 : 70-sec On - 10-s Off scans with spline smoothing 
(case 2) with the total telescope time of 400 s. CI and C2 : case-2 scans with the telescope time of 1200 s. 
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1 2 5 10 20 60 100 200 500 

Integration time [min] 

Fig. 10. SDs of the resultant spectra derived by ON-OFF scans as functions of total telescope time. 
Filled circles indicate results from sets of conventional 30-s ON + 30-s OFF scans without spectral 
smoothing. Opened squares represent those of 70-s ON + 10-s spline-smoothed OFF scans. The solid 
lines indicate median values of SDs obtained after baseline-subtraction was applied. The median was 
taken from a number of sets in 10-hour observations. The grey shade indicates maximum and mini- 
mum values of the sets. While the power-law index of the median values (solid lines) in conventional 
scans were —0.50 between 60 and 32400 s, that in proposed scan was —0.46 between 80 and 32400 s. 




o _ 
Oi 

™ — \ \ 1 1 1 \ r 

2 5 10 20 50 100 200 

Spline Node Inten/al [ch] 

Fig. 11. Standard deviations (SDs) of the OFF-source-subtracted spectra. Filled circles and opened 
squares stand for scan patterns of case 2 with the telescope times of 400 s and 1200 s, re- 
spectively. Solid and dashed lines indicate SDs before and after baseline correction with 45-ch 
spline smoothing is applied, respectively. The horizontal dotted line shows the SD level of case 
1 (conventional ON-OFF scans without spectral smoothing). Median values are plotted with 
the symbols. The grey shades indicate the peak-to-peak range of the baseline-subtracted SDs. 
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